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ABSTRACT 

Linear acceleration emission occurs when a charged particle is accelerated parallel to its velocity. We evaluate 
the spectral and angular distribution of this radiation for several special cases, including constant acceleration 
(hyperbolic motion) of finite duration. Based on these results, we find the following general properties of the 
emission from an electron in a linear accelerator that can be characterized by an electric field E acting over 
a distance L: (i) the spectrum extends to a cut-off frequency Ha> c /mc 2 ~ L(£ , /£'schw) 2 Mc> where £schw = 
1.3 x 10 18 VirT 1 is the Schwinger critical field and Ac = ft/mc = 3.86 x 1CT 13 m is the Compton wavelength 
of the electron, (ii) the total energy emitted by a particle traversing the accelerator is ^afha> c , in accordance 
with the standard Larmor formula, where <?f is the fine-structure constant, (iii) the low frequency spectrum is 
flat for hyperbolic trajectories, but in general depends on the details of the accelerator. We also show that linear 
acceleration emission complements curvature radiation in the strongly magnetized pair formation regions in 
pulsar magnetospheres. It dominates when the length L of the accelerator is less than the formation length 
ply of curvature photons, where p is the radius of curvature of the magnetic field lines, and y the Lorentz 
factor of the emitting particle. In standard static models of pair creating regions linear acceleration emission is 
negligible, but it is important in more realistic dynamical models in which the accelerating field fluctuates on a 
short length-scale. 

Subject headings: plasmas — pulsars: general — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Linear acceleration emission is the radiation produced 
when a charged particle is accelerated in a direction parallel 
to its velocity. The resulting rectilinear trajectory is a special 
case. In astrophysics, it is expected to occur in the inner parts 
of a pulsar magnetosphere, where the strong magnetic field 
enforces one-dimensional motion on the electrons, and there 
is a non-zero electric field aligned with the magnetic field. 
The calculation of linear acceleration emission is complicated 
by the fact that relatively large sections of the trajectory ra- 
diate coherently (in a sense that will be made more precise 
in Section |2j. In the celebrated case of 'hyperbolic' motion, 
in which the acceleration measured in the instantaneous rest 
frame of the particle is constant, the entire trajectory remains 
coherent and it is impossible to define a local photon emission 
rate. Linear acceleration emission is a potentially important 
component of the electromagnetic cascades that occur in pul- 
sars. These cascades are responsible for the materialization 
of the electrons and positrons that radiate in pulsar wind neb- 
ulae, and, ultimately, make th eir way to Earth as a high-energy 
cosmic-ray component (e.g.. lProfum o 2008). 

Early work on this problem in the context of pulsars treated 
the relativistic oscillator, in which the el ectric field is a lin- 
ear function of pos ition dWagonerlll969|). Linear solutions 
were c onsidered bv lMelrosd ( 19781) and lkuiipers & Volwerkl 
(119961) who emphasized linear acceleration emission masers. 
Recently, the focus has shifted to superluminal, large- 
amplitude electrostatic waves. The particle trajectory and 
resulting emission in these waves was c alculated using a 
semi-classic al formalism by iRowd d 19951) and has been re- 
analyzed by iMelrose et ail (120091) and lMelrose & Luol (120091) 
who found and discussed some significant differences with 
Rowe's results. 
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In this paper we compute the emissivity for linear acceler- 
ation in several special cases. The results enable us to esti- 
mate the emission in the general case, and we use this esti- 
mate to discuss the importance of the process in both static 
and dynamic models of the pair producing regions in pulsar 
magnetosph eres, which w e call " gap s". Our method differs 
from that of lMelrose etail ( 120091) and lMelrose & Luol d2009l) 
in several respects. In particular, by avoiding the use of the 
Airy integral approximation' we resolve a puzzle they en- 
countered concerning the applicability of Larmor's formula 
to linear acceleration. We also drop the restriction to su- 
perluminally propagating disturbances, which enables us to 
treat static gaps, as well as those containing spatially isolated 
electric field structures such as double layers. Although the 
method and the details of the emissivity we find are substan- 
tially different, the main conclusions concerning the impor- 
tance of line ar acceleration emission in pulsars are consistent 
with those of lMelrose & Luol f2009): it is unimportant for the 
parameters currently under discussion in static gap models of 
pulsars, but could be the dominant radiation mechanism in 
dynamic gap models. 

The paper is organized as follows: in section [2] we briefly 
review the formalism of classical electrodynamics concern- 
ing the radiation of test particles in prescribed vacuum fields, 
emphasizing the role of the photon formation or coherence 
length. In section [3] we apply this formalism to the well- 
known problem of hyperbolic orbits (motion with constant 
acceleration in the particle's rest frame) and derive a useful 
approximate representation of the angular and spectral distri- 
bution of radiated photons. In section|4]we use the method of 
steepest descents to compute an approximate emissivity for 
a simple wave form that can describe a stationary or mov- 
ing double layer. We also present an exact analytic solution 
for the spectrum emitted by a particle undergoing simple har- 
monic motion. In each case we demonstrate that the results 
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agree with simple estimates based on the spectrum found for 
hyperbolic orbits. Finally, we discuss in section|5]the role of 
linear acceleration and curvature emission in pulsar gap mod- 
els, comparing the relative importance of each process. Our 
conclusions are presented in section|6] 



2. EMISSIVITY 



Consider a particle that moves in vacuum. The spectral 
and angular distributio n of the energy U it radiates is (e.g., 
ISchwinger et"aT][T998h 
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where n is a unit vector in the direction of propagation of 
the emitted wave, and j(k, a>) is the Fourier transform of the 
current: 

j(k, cS)- J" d 3 * J" dt j(x, ?)exp (icot - ik ■ x) 

Xoo 
dfj8(f)exp[^(f)] . (2) 
oo 



Here, 



<p{t) -<±>t-k- x(t) , 
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x(t) and cJ3(t) are the position and velocity of the particle at 
time t. The quantity </>(?) can be identified as the phase of the 
emitted wave at the particle's position, given that <p = at 
t = 0. Thus, 



dU ofq 
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df I dr[nx j8(f)] ■ [n x j8(f + t)] 

oo U -co 



dwdQ 4-n 2 c J_ 

exp j-iW + icon ■ [x(t) - x(t + t)] /c) (4) 

For sufficiently large at, the exponential in © is sharply 
peaked around r = 0, which implies that the r-integral can be 
written as a function of the particle velocity and its derivative 
at time t, i.e., as the integral over the trajectory of a locally 
defined emissivity. For lower frequencies, however, this is 
not the case, and different parts of the trajectory contribute 
coherently to the emission. At any time f, the section of the 
trajectory that contributes coherently is that over which the 
phases <p(t) and <p(t + r) do not differ by a large amount; i.e., 
one can identify the time r co h over which <p{t) changes by 2n 
as a characteristic coherence t ime, or photon formation time 
(e.g lAkhiezer & Shul' gal 1981 . 

For a relativistic particle whose trajectory has a ra- 
dius of curvature R, the coherence time is approximately 
Max(y 2 /w, R/yc), and a local emissivity can be defined, pro- 
vided that the accelerating fields can be considered con- 
stant on this timescale. This is the normally encountered 
case of sync hrotron radiation. However, for both 'jitter' ra- 
diation (e.g., iToptygin & Fleishman! 1 19871: iMedvedevI l2000t 
iKirk & Revillel2010l) and linear acceleration emission, the co- 
herence time is not short compared to the timescale on which 
the accelerating fields vary. 

For a particle that experiences acceleration only during a 
finite time interval, t\ < t < h, the Fourier transform of the 
current can be written in terms of the trajectory within this 
interval: 



j(k,to) = iqc I 
Jti 



dt \d<p/dt) 
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as emphasized bv lSchwinger et all dl998l) . Specializing to the 
case of linear acceleration emission in vacuum, k — nu>/c and 
n-fi — fi/3 and the cosine fi of the angle between n and x — PI (3 
is a constant. Therefore 
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provided fi + 0. The expressions (|5]l and Q are the starting 
points of our analysis of the emissivity. It is worth noting that 
when Q is inserted into ([]]) the emissivity can immediately 
be integrated over frequency and angle to give 
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where = d/3/dt, ]f = mc{y,yf}) is the particle momen- 
tum 4-vector and t the proper time. The integrand on the 
right-hand side of (|7]) is just the relativistic form of Larmor's 
formula for the instantaneous power radiated. Thus, when 
integrated over a trajectory on which the acceleration is non- 
vanishing only within a finite time interval, this formula cor- 
rectly gives the total radiated power. It is, nevertheless, incor- 
rect to regard it as a truly instantaneous power, for example 
by associating a particular section of the trajectory with a part 
of the radiated power (e.g. . ISchwinger et alJ fl99&. Chap. 37). 

3. HYPERBOLIC MOTION 

The trajectory of a particle moving with constant acceler- 
ation (as measured in its instantaneous rest frame) is called 
"hyperbolic" because of its shape in the (jc, f)-plane (the accel- 
eration and velocity are assumed parallel). In classical elec- 
trodynamics in flat space, a hyperbolic trajectory results when 
a particle moves in a constant electric field which is parallel 
to its velocity, and also to the magnetic field, if any is present. 
The sometimes controversial literature on the problem of the 
radiation emitted by a particle on su ch a trajec t ory, go es back 
over a century — see, for example, iGinzburgl (119701) . In the 
last few decades, interest has arisen in the associated q uantum 
effects, the Unruh Effect and Unruh Radiation (e.g. .lUnruh 



1976 [lBell & Leinaasll987tlCrispino et al.l2008tlThirolf et al 
2009). At least in the classical limit, the conceptual problems 
associated with hyperbolic motion disappear if the particle ex- 
periences acceleration during only a finite time-interval, out- 
side of which it moves with constant velocity. We will assume 
this to be the case in the application to pulsars. 

Consider a particle of charge q and mass m moving along 
the .x-axis in an electric field E(x) = xE(x) which is constant 
and of magnitude E between x — and x — L and vanishes 
elsewhere. Assuming the velocity is also directed along x, the 
trajectory is 

x = a* [cosh (cr/a*) - 1] , ct = a* sinh (cr/a*) (8) 

forO < r < Tl = (a*/c)cosh _1 [{L + a*) /a*]. The length a* is 
the distance over which the electrostatic potential changes by 
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where Ac = fi/mc is the Compton wavelength and E = 
\q\hE/m 2 c 3 is the electric field in units of the critical field 
— for an electron or positron this is the Schwinger field 
£schw = 1.3 x 10 18 Vm _1 . The trajectory (0 is an exact 
solution of the classical equations of motion, including the 
Lorentz-Abraham-Dirac form of the radiation reaction force 
(which vanishes identically for < t < t£). 

In a gap near a pulsar surface, the magnetic field is suffi- 
ciently strong that electrons and positrons rapidly decay into 
the Landau ground state. Then, if the radius of curvature of 
the field lines is large (see section|5]) the particle motion is ap- 
proximately one-dimensional. The component along B of the 
electric field induced by rotation of the star accelerates parti- 
cles in a trajectory that can be approximated by ©, whereas 
the perpendicular component produces an E x B-drift that can 
be transformed away. Static models of gaps assume the (par- 
allel) field vanishes below the surface (x = 0) and above a 
"pair-production front" at height L, where photons produced 
in the gap create a sufficient number of charges to screen the 
field. Under pulsar conditions, one generally expects £ « 1 
and LI a* » 1 (see Sect. [5}. 

A particle that starts at rest at the surface and is accel- 
erated upwards achieves a Lorentz factor y — (x + a*)/a* 
at height x. The energy radiated in passing from the sur- 
face to the pair production front can be computed from Q 
or (O. For emission at n ■ 5t = n, a simple approximation 
can be found by first transforming into the particle rest frame 
at proper time to = (a*/c) tanfT 1 (yu), when it has reached 
the position xq = a* [cosh (cTo/a*) - 1] and a Lorentz fac- 
tor jq = (xq + a*) /a*. Denoting coordinates in this frame by 
(x', f), the trajectory is 

x = a* {cosh [c(t - r ) /a*] - 1} + x , (10) 
ct' =a*sinh[c(T- t ) /a*] (11) 

For photons of frequency to' emitted at n' - 0, i. e., normal to 
x in this frame, one finds ( ISchwinger et al.|[l998l) 
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where K\ is a modified Bessel function of the second kind. 
This approximate solution involves extending the limits of the 
integration over the trajectory to t = +oo, and is accurate 
provided the endpoints in the new (primed) frame of reference 
are far from the point at which the particle is at rest: <k 
to <K t l . Returning to the frame in which the stellar surface 
is at rest, (fl2l i describes the radiation emitted at fi = /3 — 
sjy 2 — 1 ly. Therefore, exploiting the transformations a>' = 
Deo, dU' = DdU, dQ. = £T 2 dC2' where D = y(l - 0) = 
-yl -/j 2 is the Doppler factor, 
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where 9 = cos 1 fi. The result can be checked by integrating 
over all frequencies, using the identity 
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to find 
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in agr eement with the expression given by ISchwinger et alj 
d 19981) . For small argument, K\(x) * I fx, so that the low 
frequency spectrum is fiat: 
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Using the asymptotic form of K\ for large argument, one finds 
the high frequency spectrum: 
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dtodQ 4jtc a) c 
which has a cut-off at the frequency 
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As expected, for relativistic particles, most of the energy 
is radiated close to the forward direction. The method ap- 
plies only to radiation inside a small forwardly directed cone 
of opening ang le Vl" 

- jj 2 <K 1, because otherwise xo is too 
close to the stellar surface. However, the region outside of 
this cone is not expected to contain significant power. The 
approximation also fails within a very small, forwardly di- 
rected cone, because the appropriate value of to approaches 
or exceeds Tl- This is potentially a more serious shortcoming, 
since it prevents one computing the emissivity on the beaming 
cone of a particle when it reaches its maximum Lorentz factor 
Tmax = (L + a*) /a* at the pair production front. 
In fact, an explicit expression for the emission at /j. — 

Vyinax ~ 1/Tmax can be written down by evaluating the inte- 
gral over the orbit between the limits t = -oo and t - tl- 
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where I\ is the modified Bessel fun ction of the first kind 
and L \ is the modified Struve function dAbramowitz & Stegunl 
119721) . The spectrum at different angles, using this formula 
and equation <TT~4b is shown in figure Q] At high frequencies 
id » Lc/a* 2 one finds from (f20b the asymptotic expression 
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2 t2 

q L c 
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In contrast with the emission at larger angle, this spectrum 
does not exhibit an exponential cut-off, but is a power law 
oc of 2 . This behavior is clearly related to the discontinuity 
in the electric field at x = L, In reality, the pair produc- 
tion front is not a sharp boundary, but rather a region of fi- 
nite width over which the field drops continuously to zero. 
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The width of this region is determined by the distance over 
which the electric field can induce a substantial charge sep- 
aration in the newly created pairs. Assuming these are born 
with Lorentz factor y pa ; r =s p/L, where p is the radius of cur- 
vature of the trajectory imposed by the strong magnetic field 
(see section the width of the front can be estimated to be 
roughly Ax « y pa ir<z* ~ a*p/L » a*. Thus, the emission 
produced at small beaming angles 6 ~ a* /Lis unlikely to dif- 
fer qualitatively from emission at larger ones, and the entire 
pattern can be approximated as a hollow cone, with a cut-off 
frequency that increases towards the inner edge 



dU 
dwdO 



g£ [& Kl ( sin ^)] 2 forfl */L < sin0 < 1(22) 
otherwise 



Here, we have multiplied the spectrum ( [14b by normal- 
ization constant 32 /9n such that the total energy radiated 
matches that given by Larmor's formula (0: 
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If the details of the emission on angular scales a*/L 
are unimportant, a "^-function" approximation can be used. 
While the integration over solid angle cannot be found, a sim- 
ple function that approximates the angular integrated spec- 
trum is 
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The accuracy of this approximation can be seen in Figure [2] 
where the function is plotted together with the numerically 
integrated value of (1221 for a range of L/a*. 

The maximum photon energy produced can be reasonably 
approximated as 



tlU) n 



■ - LEI a* 
= LE 2 /A C 



(25) 

As we discuss in section|5] this implies that each electron that 
traverses the gap from x = to x — L produces on average ctf 
photons of frequency a> ~ <w max , where af is the fine-structure 
constant. 

4. SPECTRA FROM MOTION IN ELECTROSTATIC WAVES 

In this section we consider two special cases for which ap- 
proximate or analytic expressions can be found. The first is 
an isolated structure containing a reversal of the the electric 
field. To facilitate the analysis, we choose a particle orbit of 
the form 



x(f) = ^tanh(fi f) 

"0 



(26) 



with |/3o| < 1, i.e., the particle starts and finishes its orbit at 
rest, while its position suffers a displacement of 2c/?o/£2o- The 
maximum particle speed is c/3q. The electric field at the posi- 
tion of the particle is 



E[x(t), t] : 



mcClo -2y8otanh(Qof) sech 2 (Q f) 



[l -$sech 4 (£2oO] 



3/2 
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Of the many functions E(x, f) that can provide such a tra- 
jectory, those representing an isolated structure that is either 
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Fig. 1. — Energy spectrum, in units of q 2 /n 2 c, radiated at diffe rent 
angles using equation (14) for sin = 0.1, 0.5, 1.0 and equation )20> for 
sin 8 = a*/L = 10~ 3 . Note that the high frequency power law at angle 
8 = a* /L is an artifact of the discontinuity in the electric field (see text). A 
color version of this figure is available in the online journal. 




10 z 

u>a*/c 

Fig. 2. — Comparison of numerical integr ation of eq. I22t over solid angle 
(solid lines) with the approximate equation 1241 (dashed lines). A difference 
of a few percent is found for oj ~ Lc/a* 2 . The spectrum is again in units of 
q 2 /n 2 c. A color version of this figure is available in the online journal. 



static or moving at constant, subluminal velocity (which is 
nowhere equal to the particle speed) are perhaps the simplest 
that are physically realistic. In the static case, the electric field 
is simply 



mc n 2/jox(i 2 ~i) 



E(x) = \ " [i-/^(i-* 2 ) 2 




jtt for |Jc| < 1 



otherwise 



(28) 



where jc = Clox/cfia, but E(x,t) must be constructed using 
a simple algorithm in the moving case. Examples of these 
structures are shown in Fig. [3] If they accelerate particles to 

relativistic velocity, i.e., Jq = 1/ ^1 - f} 2 » 1, they have 

the property that the maxima of the electric field are sepa- 
rated in space by roughly Ax ss c/yo^o and have magnitude 
£max ~ 7QincQ.o/q. Note that these are not self-consistent 
structures; they have been chosen to illustrate the properties 
of linear acceleration radiation. 
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Fig. 3. — El ectro static waves that yield an electric field at the particle posi- 
tion given by (28), with yo = 10, as a function of phase f = f2o(j8 w ? + x/c) 
for three values of /3 W (and assuming q > 0). A color version of this figure is 
available in the online journal. 



From ©, allowing the limits tip 
dU 



boo, one finds 
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This integral can be approximated using the method of steep- 
est descents. As a function of complex argument z — x + iy, 
f(z) has a suitable saddle-point at x — 0, y — arccos V/?o cos 6 
(0 < y < n/2). The result is best expressed in terms of the an- 
gle 8 between the direction of propagation of the photon and 
the .x-axis, as measured in a frame moving along this axis at 
speed Pq, in which case 



cosf? = 



cos 9 - Pq 
1 -/Jo cos 
and one arrives at the expression: 
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Thus, the spectrum rises linearly with frequency up to a cut- 
off at 



<j c = 3(l +cos0) 3/2 y^Q o /4 . 



(33) 



Apart from the low frequency spectrum, which in any case 
is not reliably given by the method of steepest descents, this 
spectrum is roughly the same as that emitted by a particle 
experiencing uniform acceleration in a field E « E max = 
yZmcQo/q over a distance L — Ax — cj (yo^oX given in ( [25) . 

The second case we consider is motion in a periodic wave. 
To find an analytic solution, we choose a particularly simple 
particle orbit 



x(t) 



cos(Qof) 



(34) 



I.e., the particle undergoes simple harmonic motion in one 
dimension. The spectrum produced by a par ticle with suc h 
a trajectory has previously been calculated by lSchottl (119121) . 
The electric field at the position of the particle in this case is 



E[x(t),t] = 



-Po cos(Q f) 



1 [[l-^ s in 2 (O 0] 
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which, as in the previous case, has maxima of the electric 
field separated in space by Ax * c/joCIq with magnitude 
Emax ~ TntncClo/q. Thus, if the particle becomes highly rel- 
ativistic, the field is large, but is confined to a small spatial 
region. Since the motion of the particle is periodic, the spec- 
trum consists of harmonics of the oscillation frequency D.q. 
The resulting power, averaged over one oscillation period, is 



2^-mQo) — 



dCldco 



where 



dP m ^Q 2 tan 2 
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(37) 



which clearly demonstrates the discrete na ture of the sp ec- 
trum. An alternative derivation to that of iSchotJ dl912l) is 
given in the appendix. Equation d37] i represents the average 
power radiated, in the m-th harmonic, into a unit solid angle 
at an angle 6 to the oscillation axis. For m < m c (see below) 
the power radiated in each harmonic at a fixed angle is propor- 
tional to m 4 ^ . Fig.|4]shows a plot of this function for different 
emission angles. 

When to is large, the distance Ho between neighboring har- 
monics becomes comparatively small, and the spectrum can 
be considered effectively as a continuum. The transformation 
to a continuous frequency spectrum can be made by formally 
replacing the summation in equation (f36t by an integral 



Z 



dm 



Taking the asymptotic value of the resulting expression, we 
find, in terms of the angle 6 measured in a frame moving at 
speed Po 



dcodQ. 



period 



sin 6 — exp 

\6ttc co c 
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where we have multiplied by a factor 2^/Qo to give the aver- 
age energy emitted in each period, The critical frequency in 
this case is 
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This result is very similar to the approximate expression for 
an isolated structure, provided the strength of the electric field 
and its spatial extent coincide. 

Approximate expressions for the emi ssion produced by a 
relativistic oscil l ator h ave been found by W agoned i 19691) and 
iMelrose & Luol (120091) . These share several properties with 
the simpler cases presented here. To compare with the emis- 
sion from a hyperbolic trajectory, note that the electric field 
varies only linearly over the wave period, so that its effec- 
tive spatial extent is Ax » c/D.q. The effective field strength, 
on the other hand is E max as y m . ix mcQ.ol q. According to 
( l25l l a cut-off is expected at w cr ; t = y 2 nax Qo, as found by 
the above authors. The low frequency spectrum, however, is 
more sensitive to the wave form, with the oscillator having 
df//d<jxlQ oc a), whereas the low frequency spectrum radiated 
by the hyperbolic orbit is flat. 

5. APPLICATION TO PULSARS 
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Fig. 4. — Envelope function for the power radiated in the m-th harmonic for 
different angles in units of q 1 Qij2nc for an electron with maximum Lorentz 
factor yo = 10, c.f. Equation {37). A color version of this figure is available 
in the online journal. 



Pulsars are thought to lose their rotational energy in a 
relativistic wind that contains a large number of electron- 
positron pairs (for a review, see lKir k et al. 2009[). These pairs 
presumably materialize in an electromagnetic cascade initi- 
ated by electrons that are accelerated in the strong induc- 
tive electric fields surrounding the star. Recent observations 
indicate that ga mma-ray pulses are unlikely to be produced 
close to the star (fAbdo & the Fermi LAT collaborationll2009t 
lAbdo et ai1l2009l) . However, this is not necessarily the region 
where the pairs materialize, the location of which is still un- 
known. Detailed models of the electromagnetic cascades hav e 
been constructed (reviewed in iHardirigl 120091 : IChengl 120091) . 
but these face difficulties in accounting for the o bservations 
dAbdo et alj|2009t Ide Jagerll2007t lTimokhinll2009h . 

The radiation processes thought to be important in pulsar 
gap models are curvature radiation and inverse Compton scat- 
tering (both resonant and non-resonant). Linear acceleration 
emission is closely related to curvature radiation, since they 
both depend on the configuration of the large-scale electro- 
magnetic fields anchored in the star. Inverse Compton scat- 
tering, on the other hand, depends not only on the local mag- 
netic field strength, but also on the target photon population. 
Consequently, to assess the possible role of linear accelera- 
tion emission, we compare its spectrum and luminosity to the 
standard expressions describing curvature radiation. 

The acceleration region or "gap" can be characterized by its 
height L, and the potential drop y max mc 2 /e across it. Assum- 
ing the electric field to be uniform, it can be expressed in units 
of the Schwinger field as 

E = A c y nwx /L . (40) 

Denoting by p the radius of curvature of the magnetic field 
lines in the gap, the characteristic or cut-off frequencies of 
linear acceleration emission, <yii nac = rnc 2 u)n nac /h, and curva- 
ture emission, w culv = mc 2 co cnrv /h, are 

maze = ymax-E = y^Xc/L , (41) 
&curv= rLuefc/P = ^Iinac (Tmax-L/p) . (42) 

The energy radiated as linear acceleration emission, U\\ m c = 
U\m 3C mc 2 and as curvature radiation t/ curv = U cmv mc 2 by a 



single electron that traverses the gap is 

#iinac=4a f y max £/3 = 4a 1 y 2 am A c /3L , (43) 
U cmv = 2a ir 4 max A c L/3p 2 = ±£>i inac (y max L/p) 2 . (44) 

These expressions show that linear acceleration emission 
and curvature radiation are complementary rather than com- 
peting processes. Each of them is an approximate evaluation 
of the radiation emitted by a charge accelerated along a curved 
trajectory, computed according to classical electrodynamics. 
In the case of linear acceleration emission, the approximation 
of large radius of curvature is used. In the case of curvature 
emission, the particle's acceleration parallel to its velocity is 
neglected. If the system is such that L > p/y mdx , i.e., if the 
formation length for curvature photons is smaller than the ac- 
celeration region, then the radiation is well described by a 
local emissivity, and this can be evaluated from the standard 
curvature formulas, because the change in Lorentz factor is 
small within the photon formation length. On the other hand, 
if L < p/y mdx , the photon formation length is roughly equal 
to the size of the system and it is impossible to define a local 
emissivity. In this case, it is inconsistent to use the expres- 
sions derived for curvature radiation, since they apply only 
when the formation length is limited by the curvature of the 
trajectory. When the formation length is instead limited by 
the size of the acceleration region, the expressions derived for 
linear acceleration emission are applicable. 

Note that these arguments do not apply to particles that are 
accelerated externally and are injected into the gap with y > 
Tmax- In this case, the effects of the longitudinal elec tric field 
are sm all and can be treated using a linear theory dMelrosd 
119781) . Strictly speaking, the arguments also do not apply if 
the particle is accelerated after leaving the gap — for example, 
if it continues to propagate along a curved field line without 
linear acceleration. However, in this case the local emissivity 
for curvature radiation can be used to describe the additional 
radiation. 

Note that the average number of photons emitted by a par- 
ticle as it traverses a photon formation length is af. This is 
true both for linear acceleration emission and curvature (and 
synchrotron) radiation. For a classical picture to be valid, the 
energy of the radiated photon must be small compared to the 
energy of the radiating particle. In the case of linear accelera- 
tion emission this implies E «: 1 . Therefore, since the energy 
radiated is ~ a{Ey nYdx mc 2 <K y mdx mc 2 , classical linear accel- 
eration emission cannot lead to a saturation of the acceleration 
process, in the sense of limiting y max . This is consistent with 
the fact that the radiation reaction force vanishes at all points 
on the trajectory except the end points. Classical curvature 
radiation, on the other hand, does saturate the acceleration, if 
the region exceeds a size L max given by 

imax = (p/Tmax) (<*f£) ■ (45) 

Gaps in pulsar m agnetospheres can be divided into two 
types : vacuum gaps dStur rock 197li lRuderman & Sutherland! 
11975b and gaps produced by space-charge limited flow 
dArons & Scharlemannlll979l) . The electric field is produced 
by the rotation of the strong surface magnetic field B S urface = 
SS crit with B clit = m 2 c 3 /eh = 4.414x 10 13 G. In vacuum gaps, 
it may be as high as B sur face times the rotation speed/c, and it is 
convenient to express it in terms of this value by introducing 
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the parameter e vac ~ 1 according to 

E = e VRC Br„ 



! e vac 5 x 10~ 6 [Pis/P] , (vacuum) 



(46) 



where r* = R^/r^, is the ratio of the stellar radius R* to the 
light-cylinder radius r^, P is the pulsar period measured in 
seconds, and we have used the standard estimate for the sur- 
face magnetic field from the rate of loss of rotational energy 



5*2.3 x 10 



-2 



(Ms) 



1/2 



(47) 



where P\s is the rate of change of pulsar period times 10 15 . 

The presence of space-charge limits the electric field 
to much smaller values. W hen the effects of frame- 
dragging are taken into acc ount dMuslimov & Tsyganlll992t 
iMuslimov & H arding 19971) one finds 



£~e space 0.1B< 



lO-WpV2p-3/2^ ( space . c harge) 



(48) 



where, in analogy with ( 1461 , we have introduced the parame- 



ter ^space ~ 1 ■ 



The ratio Ly max /p of the system size to the formation length 
for curvature photons determines which radiation mechanism 
is appropriate in a particular gap. This condition can be for- 
mulated in terms of a critical Lorentz factor: 



Tcrit = y]Ep/Ac 

{ 2.4 x 10 7 (e vac p.) 1/2 P 1/4 P\i 4 (vacuum) 

\!/2 



(49) 



\ 1.1 x 10 s (e space p,) P- 1/4 P\' 5 4 (space-charge) 

where we have written the radius of curvature of the field lines 
in units of the light cylinder radius, p = p t r^, and have in- 
serted the expressions (l46l l and (l48l for the accelerating fields. 
If Tmax > Tcrit (Tmax < Tcrit), curvature radiation (linear accel- 
eration emission) is the appropriate description. Alternatively, 
the curvature photon formation length L cmv can be evaluated: 



=p/ycrit 

f200(e vac p,) 1/2 P 3/4 P:, 1/4 cm 



(vacuum) 



4x 10 4 (e space p„) 1 P 5/4 P 15 1/4 cm (space-charge) 



(50) 



According to ( 1501 ). curvature radiation is an appropriate de- 
scription of the radiation from a vacuum gap (space-charge 
flow limited gap) provided its height exceeds roughly 200 cm 
(4 X 10 4 cm). Note, however, that these conditions are relaxed 
(in the sense that the critical gap size is reduced) if the ra- 
dius of curvature of the field lines is small compared to the 
light-cylinder radius, and also if the accelerating field fails to 
reach the full value given b y (|46|| or d48b — both of which 
were assumed in the original lkuderman & Sutherland! (1 1975b 
model. 

In current static models, the assumed gap heights greatly 
exceed L CU1V , so that linear acceleration emission does not play 
a role. Recentl y, however, time-dependent gap mode ls have 
been proposed dLevinson et al.ll2005l : lTimokhinll2009l) . These 
models assume accelerating fields that are of the order of, or 
larger than the vacuum field given in (l46b . In addition, time- 
dependent screening leads to structure on a length scale deter- 
mined essentially by the electron inertial length ( e = c/a) pe , 



where w pe is the local plasma frequency. Estimating the elec- 
tron/positron charge density as roughly equal to the charge- 
density required to screen the electric field (the "Goldreich- 
Julian" density) implies 



w pe = 1.5x lO'V-'^V 1 (51) 



so that 



4 fl0- 2 (e vac p«)- 1/2 p- 1/2 vacuum) 

— — = \ _<;/ \-V2 _i (52) 

Lcmv [5 x 10 ^e space p«j P '(space-charge) 

Therefore, linear acceleration emission is an essential ingre- 
dient in these models, because the formation length for cur- 
vature radiation photons substantially exceeds the anticipated 
length-scale of the structure in the accelerating electric field. 

6. CONCLUSIONS 

The radiation produced by a particle whose acceleration is 
parallel to its velocity, linear acceleration emission, has been 
investigated for a number of special cases. Particle motion of 
this type is expected to occur in the intense magnetic fields 
close to the surface of a pulsar. The emission produced by 
these particles is conceptually difficult due to the macroscopic 
size of the photon formation length. For the particular case of 
hyperbolic motion, corresponding to acceleration in a uniform 
electric field, the formation length is comparable to the length 
of the entire accelerating region. Linear acceleration emis- 
sion is similar to both synchrotron and curvature radiation in 
the sense that a photon is emitted with a probability at when 
a particle traverses a formation length. The essential proper- 
ties of the radiation produced by a particle on a hyperbolic 
trajectory can be understood in terms of two key parameters: 
the length L of the hyperbolic part of the trajectory, and the 
length a* over which the particle increases its energy by mc 2 . 
Some generic features of the emission mechanism in different 
electric field configurations can also be described using these 
quantities. In particular, provided L and a* are interpreted 
appropriately (see section [4] for some examples) the critical 
frequency can be expressed as dL> cr i t ~ Lc/a* 2 . 

Current models of pair production in pulsar magneto- 
spheres typically assume that high energy photons that induce 
the secondary pair cascade, are produced via either curva- 
ture radiation or inverse-Compton scattering. An alternative 
to these processes is linear acce leration emission, as investi- 
gated by iMelrose & Lud (120091) . We have demonstrated that 
this process is in fact complementary to curvature radiation, 
the transition between the two regimes depending only on the 
ratio of the system size to the formation length of curvature 
photons. For existing static polar gap models, linear accelera- 
tion is unlikely to play an important role, although this relies 
on the environmental parameters being close to the canonical 
values. The limits for both vacuum gap models and space- 
charge limited models are presented in section|5] 

Recent numerical studies provide evidence that the 
above stationary models are unstabl e to perturbations (e.g. 
iLevinson et al.ll2005llTimokhinll2009l) . Several analytic mod- 
els have been developed in which particles are accelerated in 
large amplitude electrostati c waves that propagate parallel to 
the magnetic field (see e.g. iLuo & Melrose! 120081 and refer- 
ences therein). Since the characteristic length scale of the 
electric fields in these models is on the order of the electron 
inertial length, linear acceleration emission is likely to be a 
more appropriate description of the resulting radiation. 
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APPENDIX 



Here we provide an alternative derivation to that of ISchottl d 1 9 1 21) for the spectrum radiated by a particle undergoing simple 
harmonic motion. Taking advantage of the periodic nature of the orbit, Eq. (Q} can be rewritten in terms of an infinite sum 



dU q 2 tan 2 i 



dQdw 4n 2 c 

where <p — Qo? + 2nk. Using Poisson's sum formula 

k 



i#-(^-/^Socos0) 



tc=oo ni=oo I v 

k=-co m=-oo \ w / 



and noting that u is positive, we find 



dU q 2 &l tan 2 



•n=°° ~n 

y 5(a) — raQo)ra I d<p& 



im({/)-ii/3o cos 0) 



dQdoj 4n 2 c 

The integral in this equation can be expressed in terms the Bessel functions 

dU q 2 Cl 2 tan 2 6 



6(a> - mQ.o)mi m J m (pfiom) 



m=\ 



dfidw 

where we have exploited the integral form of the Bessel function 

i m J m (t) = zt*® 1 

Finally using the identity 



6 2 (x) = lim — 6(x) 



the power radiated in the ra-fh harmonic is 



dp i du V? dP m 

= hm = > o(a> - ratio) 

dfidw t^oo T dQdw ^ dCl 



with 



dP m q 2 Q. 2 tan 2 6 



m J m Qj/3 m). 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



dQ 2nc 

The validity of this result can be checked by comparing against Larmor's formula Eq. @. This can be done straightforwardly 
using the identity 



2^ m 2 J m (mx) = 



x 2 (4 + x 2 ) 
16(1 -x 2 ) 7 / 2 



(9) 



and subsequently integrating over solid angle. Multiplying the resulting answer by the period of the waves gives the total 

energy radiated per oscillation 



i + 4 

r 



(10) 



It is easily verified that this agrees exactly with Eq. (j7]). 

The asymptotic expression for Bessel functions of large order and argument is (e.g lWatsonlll922i p. 249) 



1 \ 2(m-x) v 

Jm(x) -l/ " Kl/3 

n V 3x 



2 3/2 (m - x) 3/2 



3x 1 ' 2 



(ID 



with x < m. Thus, using the asymptotic expression for the modified Bessel functions with large argument, the spectrum at high 
harmonics is 



dP m q 2 n 2 



tan 2 6» 



dQ 2nc V2^VW^ 



raexp 



2 5 l\\ - hPq?' 2 

■ ; 

3 y/jjfio 



(12) 
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